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ABSTRACT 

Several past searches for gravitational waves from a selection of known pulsars have 
been performed with data from the science runs of the Laser Inferometer Gravitational- 
wave Observatory (LIGO) gravitational wave detectors. So far these have lead to 
no detection, but upper limits on the gravitational wave amplitudes have been set. 
Here we study our intrinsic ability to detect, and estimate the gravitational wave 
amplitude for non-accrcting pulsars. Using spin-down limits on emission as a guide 
we examine amplitudes that would be required to observe known pulsars with future 
detectors (Advanced LIGO, Advanced Virgo and the Einstein Telescope), assuming 
that they are triaxial stars emitting at precisely twice the known rotation frequency. 
Maximum allowed amplitudes depend on the stars' equation of state (e.g. a normal 
neutron star, a quark star, a hybrid star) and the theoretical mass quadrupoles that 
they can sustain. We study what range of quadrupoles, and therefore EoS, would 
be consistent with being able to detect these sources. For globular cluster pulsars, 
with spin-downs masked by accelerations within the cluster, we examine what spin- 
down values gravitational wave observations would be able to set. For all pulsars we 
also alternatively examine what internal magnetic fields they would need to sustain 
observable cllipticities. 
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1 INTRODUCTION 

Just under 2000 p ulsars are currently known 
(|Manchester et al [2OO3) and these provide enticing 
targets to search for continuous gravitational wave emis- 
sion. Searches for gravitat i onal w aves from known pulsars 
ijAbbott et allbOOsi . l2007al . l2010al ) estimate four unknown 
parameters: the gravitational wave amplitude ho , the initial 
phase relative to the electromagnetic pulse phase (po, the 
cosine of the orientation angle cos t, and the polarisation 
angle tp. As yet no signal has been seen from any pulsar, 
so observations provide no constraints on these parameters 
other than an upper limit on ho. However, with detector 
sensitivities improving over the next few years we will be 
beating spin-down based upper limits for many pulsars. 
This gives us the potential of direct gravitational wave 
detection (if gravitational waves provide a braking mech- 
anism that is comparable to the electromagnetic braking) 
therefore allowing us to place real constraints on these 
parameters. The ho parameter in particular will allow us 
to say something about the make up of the neutron staiQ 
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^ We will use 'neutron star' as a generic term for all types of 
compact star used in this paper, unless specifically stated. 



itself as it is related to the star's mass quadrupole moment, 
and therefore the equation of state (EoS). Some discussion 
of useful astroph ysics that cou ld be gained from measuring 
cost is given in Ijonej |2003), and accurate polarisation 
measurements could provide limits on, or allow the study 
of, different theories of gravity (see e.g. I Willi [200^ ) . but 
here we will concentrate on the amplitude and quadrupole 
measurement. It is important to stress that there is much 
that is unknown about neutron stars and whether they 
can produce and sustain large enough quadrupoles to be 
observable via gravitational waves, so explaining such issues 
is still an area of great interest. In this study we only discuss 
prospects for observing non-accreting pulsars. In accreting 
systems different gravitational wave emission mechanism 
may be present, and the links between electromagnetic 
observations and any g r avitat ional wave signal are less 
well known. IWatts et all (|2008l ) provide a good study of 
detection prospects for accreting neutron stars. 



1.1 Searches for known pulsars 

Since the start of science data taking from the current 
generation of interferometric gravitati onal wave detectors 
(LIGO, GEO 600, Virgo and TAMA [Abbott et all [iooial : 
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lAcernese et allioosl : lAndo et alllioosi ) searches looking for 
gravitational waves from a large selection of known pul- 
sars (millisecond and young pulsars with spi n frequencies 
greater than ^ 20 Hz) have been performed (|Abbott et all 
I2005l . l2007al . l2010al 'l. In the most recent analysis 116 pulsars 
were searched for using approximately a year and a half of 
data f rom each of the three LIGO detectors (|Abbott et all 
l2010al). which were operating at their design sensitivity 
I Abbott et alll2009ah . Unfortunately no signal was seen from 
any of these objects, but for the majority the sensitivity was 
still well above, by factors of 10 to over 100 times, their spin- 
down limits. The spin-down limit is set by assuming that the 
star's spin-down luminosity is equal to its gravitational wave 
luminosity i.e. all the rotational energy lost by the pulsar is 
due to radiation via gravitational waves from the Q22 mass 
quadrupole. This limit does require one to assume a mo- 
ment of inertia for the star, generally taken as the canonical 
value of 10''* kg m^ (or lO''^ g cm^ in cgs units) , and that its 
distance is precisely known. Fo r one object, the Crab pul - 
sar, this limit has been passed (| Abbott et alll2008al . [ioiOal ) 
although still no gravitational waves were seen, and for four 
others the upper limit obtained was within a factor of 10 
from spin-down. 

For known pulsars the parameter space to be searched 
over is comparatively small (position and phase evolu- 
tion are known) and long observation times can be used 
in a coherent way. This makes such searches more sen- 
sitive t han semi-targeted, or blind sear ches for similar 
sources (jAbbott et alll2007bl . l2008bl . I2009bl l3). although po- 
tentially will mi ss out on some in t eresti ng, but currently un- 
known, obiects. lKnispel fc AllenI (|2008l ) discusses the poten- 
tial strength of a population of gravitational wave emitting 
Galactic neutron stars. Here, we therefore concentrate our 
study on estimating the prospects for fully targeted searches. 

The paper is set out as follows: in ij2] we will assess 
the potential signal-to-noise ratios at which currently known 
pulsars could be observed with future detectors, review a de- 
tection statistic for these pulsars, and demonstrate the pa- 
rameter estimation capabilities of the standard search tech- 
nique at a variety of signal-to-noise ratios; in !}3]we review 
some estimates of the maximum quadrupole moments for 
neutron stars given a selection of EoS; in Uwe assess the 
potential signals, and associated quadrupoles, observable in 
future detectors based on spin-down limits and limits on 
our sensitivity, and how these limits can be thought of in 
terms of internal magnetic field strengths, and for globu- 
lar cluster pulsars limits on their spin-d own. Parts o f this 
work are similar i n scope to the review bv lOweril (|2006l). and 
the discussions in lAbbott et all ll2007bri and I Andersson et all 
(I2OIII) . 



2 ESTIMATING SIGNALS IN FUTURE 
DETECTORS 

The next (2"^*) generation of interferometric grav- 
itational wave dete ctors, such as Advanced LIGO 
faLIGO) (iHarrv et alll2010.). A dvanced Virgo (AdvVirgo) 
IjThe Virgo Collaborationll2009l ) , the Lar ge-scale Cryon genic 
Gravitational wave Telescope (LCGT) (|Kurodall2010l ) and 
GEO-HF (|Willkell200^ ). expect to provide order of magni- 
tude sensitivity improvements over current detectors, and 
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Figure 1. Design strain curves for Advanced LIGO (aLIGO) 
jLIGO Scientific CoUaborationl [2OO9I) (using the zero detuning, 
high laser power configuration), Advanced Virgo (AdvVirgo) (the 
baseline sensitivity used liere is tliat of 12 May 2009 found at 
|http : //wMMcasclna. Virgo ■ infn. it/advirgo/[ |, and the Einstein 
Telescope (ET) in both ET-B and ET-C configurations (taken 
from ,http : //www, et-gw. eu/et sensitivities) . 



offer the opportunity to beat spin-down limits for nearly 60 
pulsars (see 522}. A 3^'^ generation gravita tional wave de - 
tector caUed the Einstein Telescope (ET) (jPunturd [2OI0I ) 
is also under design study, and would offer another order 
of magnitude increase in sensitivity. This would bring hun- 
dreds of currently known pulsars into the range where we 
could beat spin-down limits (see i]2.3|l , and also may coincide 
with the completion of the Square Kilometre Array (SKA) 
radio telescope, which may give us a vastly larger number 
of sources to target. Estimates suggest the SKA may detect 
over half of the observable pulsars within the Galaxy giving 
~ 20 000 potential sources (jCordes et alll2004l ) with - 1000 
of them being millisecond pulsars and some of which could 
have large spin-down luminosities. Design strain curves for 
aLIGO, AdvVirgo and the ET in two different potential con- 
figurations are shown in Fig. [1] 

Assuming that the star is emitting at its spin-down 
limit, and that a fully coherent search can be performed, 
we have estimated the angle averaged signal-to-noise ratios 
(S /N) for all known pulsars for which a spin-down limit can 
be calculatecjfl given one year observation times, for joint 
aLIGO and AdvVirgo observations (ALV) and for the ET 
in two potential configurations - ET-B and ET-C (see i |2.3p . 

For the Crab pulsar we will estimate the S/N us- 
ing a limit based on the observed upper lim it, which at 
hf = 2.0 X 10"^^ (jAbbott et all l2010al ) is approxi- 

mately 7 times below the spin-down limit of 1.4 x 10~'^'* for 
a star with the canonical moment of inerti a and a d i stance 
of 2 kpc. This result also beats limits of jPalombal (|2000l ) 
based on the Crab pulsar's observed braking index. For the 
Vela pulsar and B0540-69 (J0540-6919), for which braking 



^ Many of the millisecond pulsars found so far are seen in glob- 
ular chisters, and accelerations within these clusters can mask 
the spin-down and even produce an observed spin-up, so in this 
section wc exclude these. 
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indices can also be reliably measured, we use up p er lim it es- 
timates based on those calculated bv lPalombal (I2OO0I). For 
Vela (using the Vela age of 11 000 years given in lAbadie et al 
I2OIII ) this gives limits 4 times lower than the spin-down 
Umit of 3.3 X 10"^", and for B0540-69 this gives a limit 
6 times lower than the spin-down limit of 5.2 x 10~^^. For 
these three pulsars our S/N estimates used known values of 
the pulsar inclination angle and gravitational wave polarisa- 
tion angles as t aken from fits to the pulsar wind nebulae in 
iNg fc Romanil | 20081 ) . As can be seen in Fig. [3] these three 
pulsars would have S/N greater than 5 for all future detec- 
tors if emitting at their spin-down, or observationally con- 
strained, upper limits. They are not included in the numbers 
quoted in the rest of this section. 



2.1 Assessing detection and parameter estimation 

We can estimate the S/N for all known pulsars for a given 
gravitational wave detector, or set of detectors, to assess our 
ability detect them and perform parameter estimation. The 
S/N, p, can be calculated from the square root of the inner 
product of the signal h with itself, p = ^/ {h\h), (the inner, 
or scalar, product for real time domain functions is given 
by {x\y) = {2/Sn{v)) x( t)y(t)dt) which given the time 
domain signal model of (e.g. iDupuis fc Woanll2005l ) 
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under the assumption that uT 2> 1, and where Sn{v) is the 
frequency dependent one sided noise power spectral densitjQ 
in Hz, T is the total observation time in seconds, TV is the 
number of samples used, and F+ and -Fx are the source po- 
sition a nd polarisation depende nt detector antenna patterns 
(see e.g. Ijaranowski et allll998l '). 

The S/N depends on the orientation of the pulsar, so 
upper and lower values can be set for the best and worst 
case orientations. The best case in terms of S/N is for a 
pulsar's spin axis to be along the line of sight (t = 0° or 
180°, or cost = ±1), which gives rise to circularly polarised 
radiation; the worst case is for the pulsar's spin axis to per- 
pendicular to the line of site (t — 90° or cost — 0), which 
gives rise to just linearly polarised radiation. The S/N ranges 
over a factor of ~ 3 between best and worst cases, and the 
angle averaged, or expected, S/N (averaging over a uniform 
distribution in cost, and <^o) is ~ 1.69 times below that 
for best case orientation value. For a fixed S/N the angle av- 
eraged value of ho is 1.89 times the value needed to give the 



This is related to the noise va,ricincG via — 2a^At, where At 
is the sample interval in the data time series. 



same S/N for the best case orientatiorQ (i.e. for the best case 
orientation a smaller value of ho will produce an equivalent 
S/N). In this section we will be asking the question "Given 
that the star is emitting gravitational waves at an amplitude 
X, what is the expected S/N it would have?", whereas in Ijl] 
we will be asking the question "Given that the signal has 
an S/N of X, what will the expected ho (or parameter such 
as ratio to spin-down limit or mass quadrupole that scales 
with ho) be?". When estimating the number of pulsars that 
are potentially observable using either an expected S/N cal- 
culated from a given ho, or an expected ho calculated from 
a given S/N, then the results will be slightly different. 

2.1.1 Detection statistic 

For the rest of this paper we will assume an S/N for a signal 
that we believe provides a good chance of detection. To do 
this we will use Bayesian hypothesis testing to provide a de- 
tection statistic with which to assess de te ction efficiency at 
differ ent S/N (see e.g. IClark et "ail l2007l or iPrix fc KrishnanI 
l2009h . We produced a Bayes factor in which the two com- 
peting models are that of data containing Gaussian noise 
and a known pulsar signal compared to data just containing 
Gaussian noise. Uniform priors were set for the signal pa- 
rameters, and a prior odds for the competing models of 1 was 
used. To assess this we produced Monte Carlo simulations of 
2000 signals, with fixed sky position, but randomly chosen 
values of cost, i/j and (/f>o, in different realisations of noise 
at a range of S/N. The approach is essenti ally using what 
is described as the B-statistic in IPrix fc Krishnan (,200!jh . 
which they show to be slightly more efficient, although com- 
putation ally more expensive, than the more widely u sed J-- 
statistic (|jaranowski et a]||l998l : lAbbott et aill2004bl ') . 

If we have a threshold false alarm rate due to back- 
ground of 1 per cent, then Monte Carlo simulations show 
that an S/N of 5 gives a detection probability of 95 per 
cent (see Fig. [2]) which we will use as our required value for 
confident detection in the subsequent sections. Real data is 
generally non- stationary and may contain interference, so 
this idealised detection statistic may be slightly optimistic, 
but should be close to reality, given well understood and 
cleaned data. 



2.2 Second generation detectors 

For aLIGO we assume two 4 km interferometers based at 
the Hanford site and one 4 km interferometer at the Liv- 
ingston site all with equivalent sensitivity and operatin g 
at the designed value (|LIG0 Scientific" Collaboratioiill2009l '). 



* Note that the angle averaged S/N for a given value of ho and the 
angle averaged ho for a given S/N are not interchangeable. This 
eomes about beeause for the first case the angle averaged S/N 
(p) oc ho{C), whereas in the second case the angle averaged ampli- 
tude {/io> ocp(l/C), where C = ([5(1 + cos2t)F+(tfe,V)]2 + 
[costFx(tfe,V)]')-'/' and (/(y)> = f^^;^^ p{y)f{y)dy. For an in- 
dividual pulsar, i.e. a single sky position, the dependence on po- 
larisation angle disappears for data spans a lot greater than a 
day, and as can be seen from Eqn. [2]the phase dependence is no 
longer present, so the orientation angle is the only one needing 
averaged over for which we use a uniform distribution between 
— 1 and 1 giving p(cost) = 1/2. 
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Figure 2. The top figure sliows tiie Receiver Operating Charac- 
teristic (ROC) curves for a variety of S/N (p) , with the p = 4 
curve showing the same form as that given by iPrix fc Krishnanl 
l|2009l '). The bottom figure shows the detection efficiency of the 
known pulsar search as a function of S/N given a background 
false alarm rate of 1 per cent. 

For Adv Virgo we use one 3 km interferometer at its design 
sensitivity (see Fig. [!}. Fig. Ogives the angle averaged S/N 
for 1 year of observations with all these detectors (3 aLIGO 
and 1 AdvVirgo, or ALV from now on) and shows that po- 
tentially 58 pulsars could be observed with an S/N greater 
than 5 (best case orientation gives 74 with S/N above 5 
and the worst case orientation gives 47 with S/N above 5). 
Some similar, bu t unpublished work, has been presented by 
ISantostasil (|2006D . 

2.3 Third generation detectors 

Some assessment of detectability of the currently known pul- 
sars, based on their spin-down limits, has previously been 
performed for the propo sed third generat i on de tector called 
the Einstein Telescope in lAndersson et all (|201lh (along with 
the potential for observing neutron stars through other grav- 
itational wave emission mechanisms), although this was just 
for the ET-B configuration (see below). Here, as above, we 
will calculate S/N for the Einstein Telescope. In Fig. Owe 
show S/N based on the sensiti v ity of the ET-B configura- 
tion, as described in iHild et all (|2008l ). with three interfer- 
ometers on one site (given here as th e current Virgo si te) in 
an equilateral triangle configuration (|Freise et alll2009l '). For 
this configuration of ET-B we could potentially see 312 pul- 
sars with an angle averaged S/N greater than 5 (or 405 for 
the best ceise and 223 for worst case). A 'xylophone' configu- 
ration (ET-C) consisting of different detectors optimised for 
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Figure 3. The angle averaged signal-to-noisc ratios for all cur- 
rently known non-globular cluster pulsars with measured spin- 
down values expected for ALV (upper panel) and the Einstein 
Telescope in ET-B (middle panel) and ET-C (bottom panel) con- 
figurations assuming emission at the spin-down limit (see text for 
exceptions). The if represents the value for the Crab pulsar, the 
• represents the value for B0540— 69, and the ■ represents the 
value for the Vela pulsar. The thick horizontal dashed line shows 
an S/N of 5. 
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different frequency bands, with comparable high-frequency 
sensitivity, but better lo w- frequency sen s itivity than ET-B , 
has also been proposed (jHild et alll20ld : [Preise et"aill2009h . 
In Fig. [3] we show the S/N based on the ET-C sensitiv- 
ity (again with the same layout as ET-B). For ET-C we 
could potentially see 648 pulsars with an orientation angle- 
averaged S/N greater than 5 (or 771 for the best case and 
528 for worst case). These numbers are around double those 
for the ET-B configuration due to the lower frequency sen- 
sitivity increase, suggesting that this may be the preferred 
configuration for known pulsar searches. However, we should 
be cautious as we w ill see in (and also as stated in 
lAndersson et alll201ll ) that many of these will have to be 
rather exotic, and therefore probably less realistic, stars to 
be seen at these levels. A further configuration, called ET-D 
jHild et alll201ll ). is now also available, although it is very 
comparable to ET-C for a wide range of frequencies, so we 
have not used it in these studies. 

2.4 Parameter estimation 

Our ability to constrain the physical parameters of a source 
is obviously dependent on the S/N, but the correlated phys- 
ical parameters we use to describe the signal can mean that 
rather unintuitively a higher S/N signal may not give the 
best constraints on ho- 

We have simulated a selection of pulsars that give best 
case and worst case S/N of around 5, 10, 50 and 100 for ALV 
(as shown in Fig. [3]), to assess the level of uncertainty on es- 
timates of their parameters assuming they are observable 
at that level. To do this we used the Markov chai n Monte 
Carlo technique described in lAbbott et ail (|2010al ') to pro- 
duce posterior probability density functions (pdfs) on the 
four parameters defined in [JT] Uncertainties on parameters 
can be assessed using a Fisher Information Matrix ap proach, 
for example as done bv ljaranowski fc KrolakI (j20ld l. but as 
they show, such a method will give biased error estimates 
overestimating the uncertainty for signals near cost = ±1, 
even for very strong signals with S/N approaching 1000. 

We calculate the appropriate signal strength needed to 
give S/N of 5, 10, 50 and 100 for both the best and worst case 
orientations (see Table [l| and produce pdfs on many reali- 
sations of the data to give an average pdf. These can be seen 
in Fig. [4] and the mean parameter values (not the peak in 
the posterior) and their la uncertainties obtained from the 
pdfs are given in Table [T] Examples of pulsars that if emit- 
ting at their spin-down limits might be observed with ALV 
at S/N of ~ 5, 10, 50 and 100 are respectively J1959-f 2048, 
J0737— 3039A (the millisecond pulsar in the famous double 
pulsar system), J0537— 6910 (an interesting young pulsar in 
the Large Magellanic Cloud, with a high glitch rate and 
large spin-down rate) and J1952+3252. The main results of 
Table[T] i.e. the relative uncertainties in ho for different S/N 
signals given the best and worst case orientation scenarios, 
are extrapolated to higher S/N and shown as the thick black 
lines in Fig. [8] (they are equivalent to the uncertainties on 
the quadrupole assuming the distance is precisely known). 

It can clearly be seen by looking at Eqn.[ll and from our 
posteriors in Fig. [4] that for a signal with cos (, Ri ±1 the ho 
and cos t parameter are highly correlated, which leads to the 
ho posterior being spread to higher values, and cos t extend- 
ing towards zero. In such a case the (po and tl) parameters 




Figure 5. The correlations between the ho and cost, and ipo and 
■0 parameters for a simulated signal with the best case orientation. 
The data points from the Markov chain Monte-Carlo (MCMC) 
and the probability contours derived from them are shown. The 
cross marks the parameters used to produce the signal. 



are also very highly correlated, with the posterior on ip be- 
ing flat over its range (see Fig. |4]| , and 00 only constrainable 
to within a 180° range i.e. the direction the star is rotating 
in can be found, but otherwise these parameters are ill de- 
fined to represent such a situation. These correlations can 
be seen more easily in Fig. [5] and they lead to the slightly 
paradoxical fact that when comparing equivalent signals in 
Fig. |4] we can actually place the tightest constraints on ho 
for the worst case orientation, i.e. lower S/N, signals! Due 
to these correlations we do not quote mean, or uncertainty, 
values on tp in Table [T] as they are meaningless, we also only 
quote mean values for 00 to show that the value found is 
oriented in the correct sense. If we were to cast the the un- 
certainties in i and i/) as a solid angle then at cos l — ±1 the 
associated error ellipses for different S/N would be entirely 
proportional to the error on t alone (i.e the error on psi stays 
constant). For a signal with S/N ~ 20 the best and worst 
case solid angles covering the ~ la probability contour in 
L — ip space are 0.35 sr and 6x 10~* sr respectively. The ratio 
of these areas scales very roughly linearly with S /N (within 
about a factor of 2) as 50p, up to S/N of ~ 100. 

In the case of pulsars like the Crab pulsar where X- 
ray observations and modelling of their pulsar wind nebulae 
are possible, there is additional information placing tight 
constraints on the orientation and polarisation angles. In- 
cluding this will break the degeneracies bet ween parame- 
ters, fo r example as perform in the search bv lAbbott et all 
(|2008ah . In Fig. |6] the fractional uncertainty on the ampli- 
tude estimate (averaged over many simulations) is shown 
as a function of cos t for simulated sources where in one 
case the orientation and polarisation angles are known (i.e. 
these are fixed at their known values and not searched over, 
equivalent to having very tight prior distributions on these 
parameters), and in the other case that they are unknown 
(i.e. the whole range of cos l and ip has to be searched over 
with uniform priors). It can be seen in Fig. [S]that for the 
case where all parameters are searched over the fractional 
uncertainty on the amplitude increases with cos t, with the 
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Figure 4. Each panel shows the average posterior probability density functions on the four gravitational wave signal parameters of ho, 
(po, cost and i/i for four simulated signals. The four signals in each panel (from top left to bottom right) correspond to the parameters 
given in each of the four sections (from top to bottom) of Table [T1 In an individual panel lines of the same colour represent signals with 
identical amplitude, solid lines represent the signals with the worst case orientation (cos i = 0) and dashed lines represent the signals 
with the best case orientation (cost = 1). The tapering off of the (po posteriors for the best case orientation at tt it it/2 is an artifact of 
the plotting and these parameters should really have abrupt cut-offs at these values. 



steepest rise between about cos t — 0.5-0.7. However, when 
cos L > 0.8 the fractional error slightly decreases again sug- 
gesting that a cosi = ±1 does not actually given the worst 
estimate of ho- This is actually slightly misleading as the 
reason for the downturn becomes more apparent from look- 
ing at the curves of fractional uncertainty given the true 
ho values (rather than the fractional uncertainty based on 
the mean ho value), and plots of pdfs for ho in Fig. U - 
at values of cost approaching ±1 the mean estimate of ho 
starts to increase relative to the true value, whilst the un- 
certainty starts to plateau, leading to the downturn. This 
means that estimates of ho when cos t = ±1 would still give 
values furthest from the true value, but due to the shifting 
of the whole posterior towards higher values rather than it 
widening. For the case where the cos t and t/) values are held 
fixed at their known values during the parameter estimation 
the fractional amplitude uncertainty behaves as one would 
expect and decreases as the S/N increases as cost tends to 
1. The fractional amplitude uncertainty, whether calculated 
with the mean or true ho value follow the same path, because 
the posterior stays symmetrical about the true value. 



Current detectors have systematic uncertainties in am- 
plitude and phase due to cal ibration in the reg ion of 10-15 
per cent and 5° respectively (|Abadie et alll2O10l ). For second 
generation detectors it is hoped to further reduce amplitude 
uncertainty to less than 10 per cent, but if we are lucky 
enough to see any high S/N (> 20) signals this may be the 
main source of error in the ho estimates. 

2.4-1 Dependence on distance 

The gravitational wave amplitude, and its uncertainties, can 
be directly measured. However, if we want to convert this 
into a physical quantity related to the star, such as the mass 
quadrupole moment (see Eqn.[3]) then the uncertainty on the 
distance to the pulsar will come into play. Current distance 
measurements for nearby pulsars come from parallax mea- 
surements, but for most others come from dispersion mea- 
sure observations extrap olated from a model of the Galactic 
free electron distribution lCordes fc Lazidl20o3 . For the ma- 
jority of pulsars distances generally have uncertainties at 
~ 2 per cent or more (even > 100 per cent) of the best 
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Table 1. Mean parameter estimates, and their uncertainties, at a variety of signal-to-noise ratios given best and worst case 
orientations. The true value of ho is given followed by the average (over many trials) mean estimate and standard deviation 
on each parameter. The fractional uncertainty on ho is also given. In the low S/N cases where the signal peaks near zero a 
95 per cent confidence upper limit is given. The true values of ipo and ip in all cases are tt rads and rads respectively. The 
parameter estimates have been performed assuming the ALV combination at design sensitivity. 
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Figure 6. The fractional uncertainties on amplitude (standard 
deviation over mean - cr^g/ho) from many simulations of perform- 
ing parameter estimation on a signal (with an S /N of 50 calculated 
for the worst case orientation) when all parameters are unknown 
(blue) compared to when the cos t and -0 parameters are known 
and fixed in the search. Also shown as dashed lines are the frac- 
tional uncertainties given the actual signal ho values rather than 
the mean recovered value. 



fit values. In the future, with instruments like the SKA, di- 
rect parallax measurements may push these errors down to 
the level of < 20 per cent for many millisecond pulsars out 
to 9kpc (|Smits et alll20lil ). Here we will look at how differ- 
ent distance uncertainties affect estimates of the quadrupole 
moment for a variety of S /N and for the best and worst ori- 
entations. 

It can be seen from Tableland Fig.[7]that for the low- 
est S/N signals distance errors of 10-20 per cent give only a 
relatively small increase in the uncertainty with which the 



quadrupole moment can be measured. The convergence of 
the uncertainties on the measurement for the best and worst 
case orientations, i.e. when the uncertainties become dom- 
inated by the distance uncertainty rather than that on the 
measurement, can be seen in Fig.|8] For strong signals (S/N 
> 50) the uncertainty on distance will start to dominate 
quadrupole estimates if it is > 10 per cent and will also 
swa mp any diffe rences due to orientation. 

ISetd (|2005l ) suggests that gravitational wave observa- 
tions of pulsars could be used to determine their distances to 
< 10 per cent, but only for nearby stars (< 100 pc) with very 
large S/N and high frequencies. This is therefore unlikely to 
help constrain distances better than other techniques for the 
vast majority of known pulsars. 



3 EQUATIONS OF STATE 

As discussed in lLattimer fc PrakashI (|2007h there are many 
ways to attempt to observationally constrain the EoS of a 
neutron star through inferences about their mass and ra- 
dius e.g. via various observations of binary and accreting 
systems. Observations of gravitational waves from merging 
neutron stars (or merging neutron stars and black holes) 
can also give mass and radius measurements, and informa- 
tion abo ut the EoS properties from tidal effects during the 
inspiral ([Flanagan fc Hi iidcrcr 200^) and the point at which 
the star breaks up (| Andersson "et alll201ll ) . 

However, emission of gravitational waves from indi- 
vidual neutron stars, either through a continuous emission 
from a susta ined triaxiality, or short burst s from vibrational 
modes (e.g. lAndersson fc Kokkotad Il998l ). could also pro- 
vide constraints. Due to the low amplitude of gravitational 
waves expected from these sources, such observations would 
most likely only be available for Galactic sources. Rather 
than studying explicitly what observations of continuous 
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Table 2. The uncertainties on the quadrupole moment at a variety of S/N given best and worst case 
orientations v^ith uncertainties on the distance measurement of ar = 1, 5, 10, 20 and 50 per cent. 
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gravitational waves from a triaxial neutron star can tell us 
about the stars EoS, we flip the question and look at the 
maximum magnitudes of signals we could possibly expect 
for a variety of EoS. The maximum values we discuss are 
likely to be at the extreme end of possibility and in real- 
ity they could overestimated (or possibly, but less likely, 
underestimated) by several orders of magnitude due to the 
uncertain physics of these objects. 



3.1 Quadrupole moments 

The observable in known pulsar searches that relates to 
the neutron star EoS is the gravitational wave amplitude 
ho. This is directly related to the star's I = m = 2 mass 
quadrupole moment via 



Q22 ~ ho 



\/l5/87r, 



(3) 



where r is the distance in m and u is the rotation frequency 
in Hz. The quadrupole moment is often parametrized in 
terms of the star's principal mome nt of inertia , Izz, and 
ellipticity, e, via eizz = \/87i-/15Q22 (|Owenll2005l ). However, 
these two parameters cannot be disentangled from observa- 
tions, although the range of moments of inertia can be im- 
plied from various theoretical EoS mass and radius relations 
(generally thought to be in the range 1 — 3x 10'^* kgm'^ e.g. 
lAbbott et alll2007ah , allowing ellipticities to be implied from 
observations. Often in gravitational wave literature results 
are cast in terms of ellipticity as it feels a more physically 
tangible parameter, giving a physical "size" of the distortion 

of a star. 

F o llowin g for example lUshomirskv et all (|2000l ) and 
lOwenI (|2005l ) the quadrupole for an incompressible star with 
a thin crust can be written as 



Q2 



'yfj.R^{a22) 
GM 



(4) 



where (a"22) is the weighted average strain on the crust (con- 
tributing to the 22-quadrupole), R and M are the star's ra- 
dius and mass, and 7 « 13. Often this equation is quoted as 
an upper limit on Q22 by inser ting the maximum brea king 
strain (Tmax, however as noted in lUshomirskv et all (|2000l ) the 
equation is equally applicable to providing an estimate of the 
crustal strain. The parameter in Eqn.|4]that is most depen- 
dent on the detailed make-up of the star is the shear modulus 
H, which between theories can potentially vary by many or- 
ders of magnitude. Below we will make the assumption that 
the stars are incompressible and therefore use Eqn.Uin cal- 
culations, along with an assumed 7 = 13 (for a conventional 
neutron star 7 can be related to the ratio of the thickness 
of the star's crust, AR, to it's radius via 7 « 120AR/R). It 
should also be noted that these assume the Cowling approxi- 
mation a nd neglec t self-gravity of the density perturbations. 
[Ushomirskv et all ((2000 ) showed that including self-gravity 
can increase the quadru pole calculations by between 25-200 
per cent, and similarly iHaskell et all (|2006l ) found that self 
gravity can affect results by factors of between 0.5-3. 

For the millisecond recycled pulsars spin-down argu- 
ments alone tell us that their quadrupoles must be relatively 
small (< 10^'^kgm^). Such a quadrupole would be obtainable 
with any EoS (see below) meaning that if detected they are 
not helpful differentiating between theories, although multi- 
ple detections could build up useful statistics on their prop- 
erties and limits on their internal magnetic fields. 

Here we will review some of the work presented bv lOwenI 

(|2005h regarding maximum sustainable quadrupoles for a 
variety of stellar EoS. 



3.1.1 Normal neutron stars 

For stars made from no rmal neutron star matter (neutrons, 
protons and electrons) lUshomirskv et al (|2000l) provide a 
detailed model of the quadrupole (see Eqn. 69). OwenI (|2005l ) 
applies standard numbers in this definition (and corrects the 
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Figure 7. The pdf for the quadrupole moment, Q22, given a 
simulated signal for the pulsar with ho = 7.31 xlO~^® in Table[T] 
at both cos t = (top plot) and cos t = 1 (bottom plot), with S/N 
of 50 and 141.5 respectively, over a range of distance uncertainties, 
ar, of 1, 5, 10, 20 and 50 per cent. 



definition of tlie sliear modulus to be 4 x 10^® erg cm ^, or 



2.5 X 10"" MeV fm"^) to give 



Q22 = 2.4x10^' kgm"(ao.i)i??o Ml" 



(5) 



wliere (cro.i) is an averaged strain of 0.1, Rio is tiie radius 
in units of 10km, and Mi. 4 is the mass in units of 1.4 Mq. 
Uncertainties in tlie star's mass will only affect this estimate 
by small amounts e.g. given theoretical^] a-nd observational 
bounds on the mass between ~ 1 — 2.5 Mq the quadrupole 
will only vary within about a factor of three from 0.5-1.5. 
However, for the radius, with its far larger exponent, small 



^ Theoretically the lower bou nd on neutron star mass could 
be as small as ~ 0.1 M0 jLattimer fc PrakashI l200j) . but 
we will assume our population of known pulsars is simi- 
lar to_the^ones_with obser ved masses - see e.g. Fig. 3 
of iLattimer fc PrakashI ||20071) . or the figure maintained at 
|http: //www, ste llar collapse . org/nsmasses A discussion of 
gravitational waves from neutron st ars at the lowes t end of the 
possible mass range can be found in iHorowitj \20idi ). which sug- 
gest these are potentially strong sources, but it is very unlikely 
that any known pulsar would be of this type. 




Figure 8. The uncertainties on the estimate of Q22 given in 
Table[2]for a variety of distance uncertainties and the best (dashed 
lines) and worst (solid lines) case orientations. 



differences can give a larger range of possible quadrupoles. 
If we take a theoretical range from 10-15 km then this can 
change the quadrupole by about an order of magnitude. The 
most massive stars will also have the smallest radii, so from 
this we get an uncertainty range on the quadrupole from 
the unknown mass and radius of between ~ 0.5-20 times 
the value in Eqn. [SI with the most massive, but smallest 
stars at the lower end and vice versa. 

Here we will assume the breaking s train is at the max- 
imum value of (Jmax ~ 0.1 calculated by [Horowitz fc Kadaul 
(|2009l ) (much previous work has assumed a maximum break- 
ing strain of 10~^ ^ cTmax ^ 10~^). This value of the break- 
ing strain was calculated for normal neutron star matter, but 
for other situations it may well not be a valid assumption. 
Using the higher value from the mass/radius uncertainty 
as an upper limit, and inserting in the maximum breaking 
strain, we could get normal neutron stars with quadrupoles 
of Qmax ~ 4.5 X 10^'^ kgm^. Converting this to an approxi- 
mate (order of magnitude) estimate of the ellipticity, assum- 
ing the canonical moment of inertia, would give s ~ 6xl0~^. 
Although the reasoning behind it is quite different (i.e. just 
coming from plugging in masses and radii at the extent of 
their ranges, with the maj ority of the increase over the fidu- 
cial value in iQwenI (|2005l ) coming from using a maximum 
radius of 15 km) this value is very similar to that produced 
by the perturbative a pproach to the problem performed 
by iHaskell et all (|2006l ) . Assuming the maximum breaking 
strain of 0.1 t hey would produc e a maximum quadrupole 
(see Table 4 of lHaskell et~dll2006h of Q^ax = 3.1x10^^ kg m^ 
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for a star with a mass of 1.4 M©, radius 12.3 km and crust 
thickness 1.5 km. 

3.1.2 Hybrid crystalline colour- superconducting star 

InlKriippel fc SedrakianI (|2009l ). lHaskell et all (|2007f ) andllinl 
l|2007l ) the quadrupole is calculated for crystalline colour- 
superconducting (CSS) hybrid stars. In these stars the grav- 
itational wave emission mainly comes from a deformed in- 
terior core of quark matter. The quadrupole can again be 
approximated by Eqn. but with a shear modulus given by 
l|Mannarelli et alllioOTi ') 

M = 2.47 MeVfm-^A?o/i' 400, (6) 

where Aio is the gap parameter in units of lOMeV, and 
fjLq 400 is th e quark chemica l poten tial in units of 400 MeV 
(which in iMannarelli et all |2007| is estimated to be in 
the range 350 MeV ^ ^ 500 MeV), and the stellar 
mass and radius are replac ed by those of the quark core. 
iKnippel fc SedrakianI (|2009l ) assume a range of gap param- 
eter^ 10 ^ A 50 MeV, and find a maximum core mass of 
0.8 Mq and a maximum core radius of 7 km. For the reasons 
stated in t|3.1.1l we will assume a slightly more conservative 
maximum breaking strain than that for normal neutron stars 
of (Tniax ~ 10~^ . This gives a maximum quadrupole (for A = 
50 MeV, ^lq = 500, /i = 97MeVfm"^ ^ 1.5 x 10^* J m"^') of 
Qmax ~ 1.4x10'^^ kgm^, or almost three orders of magnitude 
larger than a normal neutron star, which is not necessarily 
surprising since defor mations are in the h igh density core 
rather than the crust. iHaskell et a 1 (|2007l ) note that for a 
star with a fluid envelope around the core the quadrupole 
will be suppressed, particularly if there is not a substan- 
tial change in density when transitioning between the core 
and envelope. Converting this to an approximate estimate of 
the ellipticity, by assuming the canonical moment of inertia, 
would give an equivalently large e ~ 0.02! 

3.1.3 Hybrid and meson condensate stars 

lOwenI (|2005l ) also looks at hybrid stars with charged meson 
condensates and quark-baryon cores. For these the shear 
modulus is given by 

/i = 0.25 MeVfm-^gL-D?5S3"o^ (7) 

where go. 4 is the charge density of quark droplets in units of 
-0.4e, Di5 is their diameter in units of 15 fm, and S30 is their 
spacing in u nits of 3 fm. Following the correction for charge 
screening m lOwenI (|2005l ) an upper limit on shear modulus 
IS given as fi ~ 1.3xlO-^MeVfm'^ If we evaluate Eqn. [4] 
with this, again substituting for the core radius (which as in 
tj3.1.2l we will set as 8 km) and using a fiducial stellar mass 
of 1.4 M0 and a maximum breaking strain of (Jmax ~ 10"^, 
we get an upper limit on the quadrupole moment of Qmax ~ 
3.5 X 10''^ kgm^. This is less than t he extremal value for a 
normal neutron star. The work bv iKurkela et U (l2010al lbh 
into cold quark matter suggest that hybrid stars (with pure 

° iMannarelli et a 3 ||2007|'I gives a ran ge from 5 A 25 MeV, 
although iKnippel fc SedrakianI \200^ ) suggest that for the low 
temperature CCS phase appropriate for these cores larger gap 
values, maybe up to 100 MeV, are possible. 



phases of hadronic and quark matter) could have masses up 
to --^ 2.1 Mq and radii of ~ 13 km although these are likely 
to have only a tiny quark core i.e. they will mainly look like 
a normal hadronic neutron star. However, they show that 
stars with mixed phases of quarks and hadrons can have 
masses of up to ~ 1.9 Mq and radii of ~ 11 km, which with 
the above assumptions allows larger quadrupole moments 
of Qmax ~ 1.8 X 10^^ kg m'^, comparable to the maximum 
obtainable for a normal neutron star (although requiring a 
smaller breaJting strain). Converting this to an approximate 
ellipticity, assuming the canonical moment of inertia, would 
give e « 2x10"''. 



3.1.4 Solid strange stars 

In lXul (|2003l ) the possibility is presented th at neutron stars 
could be made of solid strange quark matter. [Xul (|2003l ) uses 
the observation of kHz quasi-periodic oscillations (QPOs) in 
X-ray bursts from neutron stars in X-ray binaries, and their 
potential association with torsional modes in the star, to give 
a shear modulus for solid s t range stars of ~ 4x10"^^ J m^'^. 
As pointed out by lOwenI (|2005l ) the identification of the 
QPO frequencies with torsional modes is somewhat prob- 
lematic, but similarly we will take this shear modu lus to 
be an upper limit for strange stars. Also iLirJ (|2007h notes 
that the theoretical arguments for solid strange stars are 
less robust that those for the crystalline colour supercon- 
ducting stars discussed in H3.1.2I The masses and radii for 
models of st range quark stars ca n be s een, for example, 
in Fig. 2 of iLattimer fc PrakashI ^0^. Although these 
may not necessarily hold for solid strange stars we will 
use the range of masses from 0.5-2 Mq and radii from 8- 
11km given by this figure to estimate the likely range of 
quadrupoles from Eqn. 2) Unlike a norma? neutron star mod- 
els of quark stars show a decrease in mass with a decrease in 
radius. Using the maximum mass and radius and a break- 
ing strain of cr = lO"'^, gives the highest quadrupole of 
Qmax ~ 3.5 X 10'''' kgm'^. Converting this to an approximate 
ellipticity, assuming the canonical moment of inertia, would 
give £ 5x lO""*. 



4 LIMITS FROM DETECTION 

The spin-down limit assumes that the total spin-down lu- 
minosity of the pulsar is emitted as gravitational radiation, 
but if we assume that the currently known pulsars are only 
just observable at an S/N of 5 we can estimate what angle 
averaged percentage of the spin-down power is going into the 
gravitational wave emission. This is shown, for ALV, ET-B 
and ET-C in Fig.[9| We can see that for the majority of pul- 
sars for which we could beat the spin-down limit more than 
1 per cent of the spin-down power would need to be lost via 
gravitational waves for us to just observe them (with only 12 
per cent, 26 per cent and 34 per cent requiring less than this 
for ALV, ET-B and ET-C respectively). Less than 1 per cent 
would be consistent with the one pulsar for which we have a 
reliable limit on this for (the Crab puls ar) that is losing les s 
than 2 per cent of its power in this way (j Abbott et al|[2010al ). 
However, for the majority of pulsars we have no reliable way 
to estimate how much energy is lost via gravitational waves. 
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Figure 9. The percentage of the spin-down power of Icnown pul- 
sars emitted via gravitational waves if the pulsar is just observable 
with an S/N of 5. Squares represent observation with ALV, circles 
represent ET-B and stars represent ET-C. 



especially given that only a few young pulsars (and no mil- 
lisecond pulsars) have measured braking indices, so 10s of 
per cent are not ruled out. 

As we saw in ^Hthere are large uncertainties in many of 
the parameters that give rise to a neutron stars quadrupole 
moment. The stars could also just be intrinsically quite un- 
strained i.e. {(J22) may just be small. This, and the fact that 
many parameters are highly correlated, means that gravita- 
tional wave observations from the quadrupole of a triaxial 
neutron star will generally be unable to pin-down much of 
the physics giving rise to it. However, as we have seen in 
^ there are reasonably hard upper limits on the emission 
allowable for stars with different EoS, so if we see emission 
above a particular EoS limit then these observations could 
rule that out for a particular pulsar, or alternatively rule 
others in. 

With this in mind we will look at all currently known 
pulsars within the expected sensitive bands of ALV and the 
Einstein Telescope (in ET-B and ET-C configurations) and 
say whether they could be constrained to a particular EoS, 
or whether a variety of EoS are valid. We will do this for 
pulsars that the spin-down limit suggests are observable with 
an S/N of 5 or greater. However, it should be noted that our 
list of EoS is not exhaustive, and other possibilities may exist 
(whether they have already been theorised or not), so here 
constraining a pulsar to a particular EoS may not be the 
end of the story. 

The quadrupole moments for each pulsar that could be 
observed with ALV, ET-B and ET-C are shown in Fig. [TO] 
both for the case that the pulsar is emitting at its spin- 
down limit and also an angle averaged scaled value that 
assumes the pulsar is only just observable at an S/N of 5 
(the scaling factor for each pulsar can be obtained by taking 
the square root of the values from Fig. [§)). Maximum values 
for the different EoS taken from ^ are also shown. The 
values are calculated from known spin-down limits, which 
are not available for globular cluster (GC) pulsars, so for 
these we have limits based on two assumptions. The first 
is a spin-down-based limit that assumes (see ^4.1^ that the 



star has a maximal gravitational wave spin-down of = 
—5x10"^'' Hzs~^, and the second is, as above, assuming that 
the star is observable at an S/N of 5. As shown in tj4. li the 
values of ho needed to produce the levels of emission required 
for this S/N can be converted into equivalent gravitational 
wave spin-down values. For ALV about half, and for ET- 
B/C all bar one of the pulsars (see Fig. [12}, are well within 
a range that could be masked by intra-cluster accelerations 
(we exclude the others). The number of pulsars with angle 
averaged quadrupoles below the spin-down limits for ALV, 
ET-B and ET-C are 59, 295 and 624 respectively (or 77, 408 
and 774 respectively for the best case orientation, and 50, 
408 and 531 for the worst case orientation). Note that due 
to the reasons discussed in i)2.1l these numbers are not quite 
the same as those suggested from an angle averaged S/N. 

From Fig. 1101 (summarised in histogram form in Fig. lll|l 
we see that for the vast majority of young pulsars (gen- 
erally those below ~ 40 Hz) the S/N limits correspond to 
large quadrupoles only supportable by the most extreme 
EoS mainly through the the application of a large break- 
ing strain. We unfortunately cannot expect the majority of 
these quadrupoles to be realistic unless there is a very unex- 
pected large population of highly deformed stars. However, 
this could be realistic for a small subset of these pulsars, or 
some as yet undiscovered sources. 

It should be stressed that these pulsars could however 
all be very smooth. The spin-down limits alone in Fig. [10] 
show that the lowest quadrupoles are only a few 10^^ kg m'^ 
(or converting to ellipticities a few 10~^) and all pulsars 
could be at or below this level. However, even at this low 
level Fig. 1101 shows that these may be detectable at S/N of 
a few with ET. 

A summary of the total numbers of pulsars below each 
EoS limit for pulsars within, and outside, GCs is given in 
Table [1 



4.1 Spin-down limits for Globular Cluster pulsars 

Several globular clusters have been specifically targeted in 
radio searches for millisecond pulsars due to having a large 
population of old stars. This has lead to a selection effect 
giving large numbers of the currently known pulsars residing 
in GCs. Pulsars within a GC will undergo significant accel- 
erations, which contaminate their true spin-down rate Vint 
via 

ail 

!>obs = Vint 1^ (8) 

C 

where i>ohs is the observed spin-down, u is the pulsar fre- 
quency and a|| is the pulsar's acceleration along the line- 
of-site - in fact many are observed to have a spin-up due 
to the accelerations (as they are millisecond pulsars they 
have small intrinsic spin-downs anyway, so it is quite easy 
for the accelerations to swamp this value). We can therefore 
use gravitational wave observations of this set of pulsars, 
either through upper limits or detections, to set limits on 
the gravitational wave component of the intrinsic spin-down 
rate via 

2 /2cV/.\ 
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Figure 10. Necessary quadrupole moments of known pulsars, 
based on spin-down limits, for them to be observable with (from 
top to bottom) ALV, ET-B and ET-C. Limits based on the ob- 
served spin-down, or a spin-down of i> = — 5 X 10~^^ Hz s~^ for 
GC pulsars, are given as crosses (blue for non-GC pulsars and 
black for GC pulsars). Angle averaged quadrupole values based 
on emission at an S/N of 5 are given as triangles. The approx- 
imate maximum quadrupoles from the various EoS discussed in 
^are also plotted. 




31 32 33 34 „ 35 



38 



Figure 11. Quadrupole moments of known pulsars fPig. I10(l. for 
them to be observable with ALV, ET-B and ET-C if emitting at 
the spin-down limit (top) and if scaled as to be just observable 
at an S/N of 5 (bottom). The filled histograms show just those 
pulsars outside of globular clusters and the unfilled histograms 
also include pulsars within globular clusters. 



We cannot however limit the overall spin-down as, unless 
the stars acceleration can be independently assessed and 
the observed spin-down calculated, there is no way to know 
what fraction of the total spin -down i s due to gravitational 
wave emission. As stated in lOwenl (|2006h by looking at 
the observed spin-downs given for GC pulsars in the Aus- 
tralia Telescope Nation al Facility (ATNF) Pulsar Catalogue 
^Manchester et ciill2005l ) it would be hard for cluster dynam- 



ics to mask spin-downs larger than v 



-5x 10"" Hzs- 



(this being the largest observed spin-up seen for any GC 
pulsar), so if gravitational wave observations can limit val- 
ues to smaller than this, then results could be providing new 
information. 

In Fig.[T2]the values of the gravitational wave spin-down 
that the GC pulsars would require to be observed at S/N of 
5 are given. For ALV 52 of the 103 observable GC pulsars, 
and for ET-B and ET-C 1 of the 107 observable pulsars, 
would need spin-downs greater than — 5x 10^^"^ Hzs^^, i.e. 
values that probably could not be masked by accelerations 
and therefore would already be seen as having large spin- 
downs from radio observations. For the rest of the pulsars 
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Figure 12. Spin-downs that would be needed in currently known 
GC pulsars for them to be observed with (from top to bottom) 
ALV, ET-B and ET-C at an S/N of 5. 



the spin-downs they would need could relatively easily be 
masked by accelerations, so they are definitely worthwhile 
as targets for gravitational wave searches. 

These limits are the smallest values that the spin-down 
could have to be observable, but for many of these pulsars 
there is a 1-2 orders of magnitude range between these limits 
and the (observationally motivated, but slightly arbitrary) 



maximum limit of ~ — 5x 10~^ Hzs~^. So, there is quite 
a lot of leeway for these pulsars to be emitting at greater 
levels than that producing an S/N of 5. Also, it can be seen 
from Fig. [10] that for the majority of these pulsars, both 
the maximum and minimum observable quadrupoles are in 
ranges allowable by all equations of state. 



4.2 Limits on magnetic fields 

Deformations of a neutron star can be supported purely by 
the strain that the star can sustairQ, although the mecha- 
nism giving rise to this strain may be unknown. However, 
magnetic fields give both a mechanism for producing strains, 
and a way of sustaining that strain. The external mag- 
netic fields of most pulsars can be estimated by assuming 
that all spin-down is due to magnetic dipole radiation, and 
these show that for millisecond pulsars the magnetic field 
strengths are comparatively small (10* to lO' gauss) and 
nowhere near enough to sustain an appreciable deformation 
on the star. The young pulsars have external dipole fields 
about 1000 times larger at ~ 10^^ gauss, but this is gener- 
ally still not enough to give deformations that would produce 
observable gravitational waves. However, these stars could 
potentially have internal fields far larger than the external 
dipole, which would be large enough to give ri se to gravita- 
tional w ave producing distortions. As done in I Abbott et all 
(|2010al ) we can use gravitational wave observations, or up- 
per limits, to place limits on th is internal fiel d strength for 
all the currently known pulsars. ICutleJ (|2002l ) predicts that 
in normal neutron stars toroidal magnetic fields could give 
rise to prolate stars with ellipticities of order 



£ = 1.6x10" 



(Bi5)G, 
(5?5>/G^ 



B < 10'" G, 
B > lO'^^G, 



(10) 



where (,6 15) is the volume av eraged magnetic field in units 
of 10'^ G. lHaskell et all (|2008l ') also study the role of internal 
magnetic fields, both entirely poloidal and toroidal, and how 
this would effect the star's ellipticity (in particular for a star 
with an EoS described by an n = 1 polytrope). They give 



i?^oA/^4^(B?.)x^2xlO-^ poloidal, 
'■^^ ' ' -lxlO-^ toroidal. 



(11) 



which for the toroidal case is similar to that o f lCutleil(|2002D . 
Very similar limi ts for toroidal fields in sup erconducting 
stars are given bv lAkgiin fc WassermanI (|2008l V These equa- 
tions can be re-arranged to give limits on the magnetic fields 
given gravitational wave observations of the quadrupole (see 
^ and an assumed moment of inertia, which here we will 
take as lO'** kg m^ (as noted above this is probably a lower 
limit and could differ by up to a factor of ~ 3). Potential 
measurements of the field strength for both the poloidal and 
toroidal cases if signals were observed at an S /N of 5 for the 
ALV, ET-B and ET-C set ups are given in Fig. [H 



^ How long such a deformation could be sustained due to visco- 
elastic creep s moothing it out is somethin g that needs fur- 
ther study, e.g. IChugunov fc Horowit j 1 I2OIOI ) who suggest short 
timescales of a few years for hot stars, but far longer for cooler 
stars (although they note extrapolations to lower temperatures 
and longer timescales must be treated with care). 
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Table 3. Summary of the quadrupole moment limits given in Figs.FlOland llll For comparison the total number 
of non-GC pulsars with angle averaged quadrupolcs below the spin-down limit and observable with an S/N 
> 5 are 59, 295 and 624 for ALV, ET-B and ET-C respectively. The total number of GC pulsars potentially 
observable with S/N > 5 and a spin-down limit smaller than — 5 X 10~^^ Hz s~^ is 51 and 106 for ALV, and 
ET-B/C respectively. 



Number of pulsars below limit 
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3 
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spin-down limit 
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90 
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53 
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99 
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97 
101 


ET-C 


spin-down limit 
S/N=5 

GC spin-down limit (i> = 
GC S/N=5 


-5x10- 


13HZS-1) 
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618 
106 
106 


131 
421 
101 
104 


43 
152 
99 
102 


41 
93 
97 
101 



We can see that for the young pulsars (with frequen- 
cies < 40 Hz) internal magnetic fields of greater than 10^* G 
would be required to observe them, which is of a similar 
strength to the external fields of a magnetars at 10^* — 
10^^ G. For the millisecond pulsars, in particular the non-GC 
pulsars, internal poloidal fields of a few lO^'^ G could provide 
observable signals, although if there are only toroidal fields 
then a couple of orders of magnitude higher would be nec- 
essary. Of course the field geometry could be complex and 
a combination of toroidal and poloidal components, and the 
calculations above rely on a specific EoS. 



5 DISCUSSION 

We have looked at our ability to detect and estimate param- 
eters from gravitational wave observations for known (non- 
accreting) pulsars. Using Bayesian hypothesis testing, we es- 
timate that signals with an S/N of 5 could be detected with 
95 per cent efficiency. This assumes that the data is Gaus- 
sian, and experience of real detector data shows that this 
assumption can be reasonable on short time scales and with 
effectively cleaned data, although there will still probably 
be some small degrading of detection ability. 

Once detected, we have shown how estimates of vari- 
ous parameters will be affected as signals increase in S/N, 
and also how the pulsar's orientation affects this. We see 
that for strong signals in the case where the orientation is 
most favourable and gives a larger S/N (i.e. the gravita- 
tional waves are circularly polarised) the uncertainty on the 
gravitational wave amplitude will actually be larger than an 
equivalent amplitude source with a worse orientation. So, 
when detections become regular the best parameter con- 
straints will actually be made for the linearly polarised 
sources. However, we also see that unless the distance to 
these pulsars can be measured to better than ~10 per cent 
then the fractional error on the quadrupole calculated from 
the gravitational wave amplitude will be relatively insensi- 



tive to the orientation, because the uncertainty is dominated 
by the distance uncertainty. 

We have seen the sorts of mass quadrupoles that would 
be necessary to observe the set of currently known pul- 
sars with future gravitational wave observatories, and com- 
pared these to maximum theoretical predictions for a va- 
riety of EoS. As has pr eviously been noted many times 
(e.g. lAbbott et al l2007al l for the majority of millisecond 
pulsars, if we are able to beat their spin-down limits and 
observe them, the quadrupoles (or ellipticities) they would 
have are sustainable by all EoS and would therefore not 
be able to constrain the type of matter in the star. How- 
ever, with observations of many pulsars useful population 
statistics could be obtained and potential differences be- 
tween different populations explored (e.g. GC pulsar and 
non-GC pulsars). To be observed with ALV, ET-B or ET- 
C the majority of young pulsars (with frequencies < 40 Hz) 
would have to have quadrupoles greater than « lO^'^ kgm'^, 
which would only be sustainable for highly strained crys- 
talline colour-superconducting stars or solid strange quark 
stars. Therefore if such stars were observed it would provide 
a great deal of insight into neutron star matter. The number 
of pulsars potentially observable at an S/N of 5, and consis- 
tent with being a normal neutron star is 11 for ALV, 103 for 
ET-B and 152 for ET-C. It is important to note that this 
does not mean that they will have these quadrupoles and 
they c ould all be smoother th an our limits are able to set. 

In lAndersson et"ail ()201lh . which studies the prospects 
of observing gravitational waves from neutron stars with 
ET (in particular the ET-B configuration), they take the 
assumption that realistic stars would optimistically have el- 
lipticities < 10~^ and conclude that only some of the 100s 
of millisecond pulsars currently known would be potentially 
observable. From this work we can look at the number of 
known pulsars we might expect to be observable at an S/N 
greater than 5 given that they are losing less than a spe- 
cific percentage of their spin-down power through gravita- 
tional wave emission, and have quadrupoles below specific 
levels (this is shown in Table [J). We find that for ALV if we 
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Figure 13. Potential internal magnetic field strengths that would 
be needed in currently known pulsars to observed them with an 
S/N of 5 with (from top to bottom) ALV, ET-B and ET-C. 



take pulsars that would only need to be emitting less than 
50 per cent of their spin-down luminosity via gravitational 
waves to be observable at S/N of greater than 5, and had 
quadrupoles less than Q22 < lO^^-lO^^ kgm'^ (or ellipticities 
of < 10~*-10~^) , then only three of the currently known pul- 
sars would be observable. We find that none are observable 
if requiring these quadrupole levels and also emitting less 



Table 4. The number of pulsars potentially observable given 
gravitational wave emission consuming less than various percent- 
ages of the total spin-down power, and requiring quadrupoles less 
than the given values. 
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than 10 per cent of the spin-down power via gravitational 
waves. For ET (in either ET-B or ET-C configurations) if 
we take pulsars which would only need to be emitting less 
than 10 per cent of their spin-down luminosity, and with 
quadrupoles < 10^°-10^^ kgm'^, then 6-11 currently known 
pulsars would be observable respectively. For the one pul- 
sar that we can currently calculate a limit on the spin-down 
power we see that less than 2 per cent of the spin-down lu- 
minosity is emitted via gravitational waves, so using this as 
a guide and taking only pulsars that could emit less than 
1 per cent of their spin-down power in gravitational waves 
and still be observable (and with the above quadrupoles) 
we find only one pulsar. These numbers are based on cur- 
rently known pulsars and use only the non-GC pulsars; for 
the millisecond pulsars that these results represent (due to 
the small quadrupole requirements) approximately half are 
within GCs, so assuming there is no difference in spin-downs 
between GC and non-GC populations of pulsars then the 
number of observable pulsars could double. Also, as stated 
in i}2l future radio telescopes such as the SKA could observe 
~ 1000 millisecond pulsars increasing the current number 
by a factor of 5. This could give 10s of observable pulsars 
for ET. 

We see that for many millisecond pulsars, if they have 
internal poloidal magnetic fields similar in strength to the 
external fields of young pulsars, or toroidal field similar in 
strength to the external field of magnetars, then they may 
sustain ellipticities that make them observable. In reality 
the internal field will probably consist of both poloidal and 
toroidal components. Young pulsar would require far higher 
magnetic fields for them to sustain ellipticities that would 
allow them to be observable. 

The above estimates have all assumed a year of obser- 
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vation. However, the future detectors could run over several 
years, and the third generation detectors may even be the 
premier gravitational wave observatories spanning decades. 
There could also be other similarly sensitive detectors added 
to the network (e.g. the LCGT). All of these would increase 
the ability to make observations of, and hopefully refine pa- 
rameter estimation for, known pulsars, although the scaling 
of these improvements (for equivalent detectors) will only 
be the square root of the observation time. Better prospects 
could come from a far larger selection of pulsar targets from 
future radio/X-ray/7-ray surveys. 

One of the main points to note here is that all the above 
calculations have significant uncertainties to them. The the- 
oretical study of neutron star EoS is still an area of study 
with large uncertainties. Neutron star mass measurements 
currently seem to favour fairly conventional EoS, and our 
upper limits for these may be a fair representation of reality. 
However, the current crop of known pulsars may harbour a 
variety of stellar types and our observations may be the only 
way to differentiate between them. More theoretical under- 
standing through simulations of maximum breaking strains 
of quark matter and hybrid star matter would be of great 
help in making these constraints. In addition to determin- 
ing whether a star is physically capable of supporting an 
observably large strain it is of great importance for future 
studies to find mechanisms of producing and sustaining such 
strains. 

The study here has also assumed emission from a triax- 
ial neutron star emitting at precisely twice the observed elec- 
tromagnetic frequency via the Q22 mass quadrupole. Poten- 
tially more interesting physics could be extracted if the grav- 
itational wave signal and electromagnetic signal were not so 
closely aligned (e.g. if pulsar timing noise is only present 
in th e electromagnetic signal I Jones! 120041 : IPitkin fc WoanI 
I2OO7I ). Deviations between the two could uncover informa- 
tion about the coupling of the gravitational wave producing 
component and electromagnetic-produci ng compone nt, and 
emission may occur at the rotation rate (|jonesll2010l ). Also, 
gravitational waves from other vibrational modes, such as 
long-lived r-modes (| Arras et all l2003: Owen 2010,;), or funda- 
mental modes (possibly excited by glitches), may be able to 
tell us far more about the neutron star EoS. 



ACKNOWLEDGEMENTS 

This work has been funded under a UK Science and Tech- 
nology Facilities Council rolling grant. I am very grateful 
to many members of the LIGO Scientific Collaboration and 
Virgo Collaboration continuous wave search group for useful 
discussions and suggestions on the issues in this paper, and 
for reading through early drafts. In particular I would like to 
thank Graham Woan, Ian Jones, Ben Owen and Keith Riles 
for their comments. I would also like to thank Roy Smits for 
the useful information he provided. 

REFERENCES 

Abadie J. et al, 2010, Nuclear Instruments Methods Phys. 
Res. A, 624, 223 



Abadie J. et al, 2011, preprint (arXiv:1104.2712v2 [astro- 
ph.HE]) 

Abbott B. P. et al, 2004, Nuclear Instruments and Methods 

in Physics Research A, 517, 154-179 
Abbott B. P. et al, 2004, Phys. Rev. D, 69, 082004 
Abbott B. P. et al, 2005, Phys. Rev. Lett., 94, 181103 
Abbott B. P. et al, 2007, Phys. Rev. D, 76, 042001 
Abbott B. P. et al, 2007, Phys. Rev. D, 76, 082001 
Abbott B. P. et al, 2008, ApJ, 683, L45 
Abbott B. P. et al, 2008, Phys. Rev. D, 77, 022001 
Abbott B. P. et al, 2009, Reports on Progress in Physics, 

72, 076901 

Abbott B. P. et al, 2009, Phys. Rev. Lett., 102, 111102 
Abbott B. P. et al, 2009, Phys. Rev. D, 79, 022001 
Abbott B. P. et al, 2010, ApJ, 713, 671 
Acernese F. et al, 2008, Classical and Quantum Gravity, 
25, 184001 

Akgiin T., Wasserman I., 2008, MNRAS, 383, 1551 
Andersson N., Kokkotas K. D., 1998, MNRAS, 299, 1059 
Andersson N, Ferrari V, Jones D. I, Kokkotas K. D., Krish- 

nan B., Read J., RezzoUa L., Zink B., 2011, Gen. Relativ. 

Gravitation, 43, 409 
Ando M. The TAMA Collaboration, 2005, Classical and 

Quantum Gravity, 22, 18, S881 
Arras P., Flanagan E. E, Morsink S. M., Schenk A. K., 

Teukolsky S. A., Wasserman I., 2003, ApJ, 591, 1129 
Chugunov A. I., Horowitz C. J., 2011, MNRAS, 407, L54- 

L58 

Clark J., Heng I. S., Pitkin M., Woan G., 2007, Phys. Rev. 
D, 76, 043003 

Cordes J. M., Kramer M., Lazio T. J. W., Stappers B. 

W., Backer D. C, Johnston S., 2004, New Astronomy 

Reviews, 48, 1413 
Cordes J. M., Lazio T. J. W., 2002, preprint (arXiv:astro- 

ph/0207156) 
Cutler C, 2002, Phys. Rev. D, 66, 084025 
Dupuis R. J., Woan G., 2005, Phys. Rev. D, 72, 102002 
Flanagan E. E., Hinderer T., 2008, Phys. Rev. D, 77, 

021502 

Freise A., Hild S., Somiya K., Strain K., Vicere A., Barsug- 
ha M., Chelkowski S., 2009, preprint (arXiv;0908.0353vl 

[gr-qc]) 

Harry G., The LIGO Scientific Collaboration, 2010, Clas- 
sical Quantum Gravity, 27, 084006 

Haskell B., Jones D. I., Andersson N., 2006, MNRAS, 373, 
1423 

Haskell B., Andersson N., Jones D. I., Samuelsson L., 2007, 
Phys. Rev. Lett., 99, 231101 

Haskell B., Samuelsson L., Glampedakis K., Andersson N., 
2008, MNRAS, 385, 531 

Hild S., Chelkowski S., Freise A., 2008, preprint 
(arXiv:0810.0604v2 [gr-qc]) 

Hild S., Chelkowski S., Freise A., Franc J., Morgado N., 
Flaminio R., DeSalvo R., 2010, Classical Quantum Grav- 
ity, 27, 015003 

Hild S. et al, 2011, Classical Quantum Gravity, 28, 094013 
Horowitz C. J., Kadau K., 2009, Phys. Rev. Lett., 102, 19, 
191102 

Horowitz C. J., 2010, Phys. Rev. D, 81, 103001 
Jaranowski P., Krolak A., 2010, Classical Quantum Grav- 
ity, 27, 194015 



Prospects of observing continuous gravitational waves from known pulsars 17 



Jaranowski P., Krolak A., Schutz, B. P., 1998, Phys. Rev. 

D, 58, 063001 
Jones D. I., 2004, Phys. Rev. D, 70, 042002 
Jones D. I., 2007, Ap&SS, 308, 125 
Jones D. I., 2010, MNRAS, 402, 2505 

Knippel B., Sedrakian A., 2009, Phys. Rev. D, 79, 083007 
Knispel B., Allen B., 2008, Phys. Rev. D, 78, 044031 
Kurkela A., Romatschke P., Vuorinen A., 2010, Phys. Rev. 
D, 81, 105021 

Kurkela A., Romatschke P., Vuorinen A., Wu B., 2010, 

preprint (arXiv:1006.4062) 
Kuroda K., The LCGT CoUaboration, 2010, Classical 

Quantum Gravity, 27, 084004 
Lattimer J. M., Prakash M., 2001, ApJ, 550, 426 
Lattimer J. M., Prakash M., 2007, Physics Reports, 442, 

109 

LIGO Scientific Collaboration, 2009, "Advanced LIGO 
anticipated sensitivity curves", LIGO-T0900288-v2, 
|https : //dec ■ ligo . org/cgi-bin/DocDB/ShowDocumeiit?docid=2974| 

Lin L-M., 2007, Phys. Rev. D, 76, 081502 

Manchester R. N., Hobbs G. B., Teoh 
A., Hobbs M., 2005, AJ, 129, 1993 
|http : //www ■ atnf . csiro . aii/researcli/pulsar/psrcat/| 

Mannarelli M., Rajagopal K., Sharma R., 2007, Phys. Rev. 
D, 76, 074026 

Ng C.-Y., Romani R. W., 2008, ApJ, 673, 411 

Owen B. J., 2005, Phys. Rev. Lett., 95, 21, 211101 

Owen B. J., 2006, Classical Quantum Gravity, 23, SI 

Owen B. J., 2010, Phys. Rev. D, 82, 104002 

Palomba C, 2000, A&A, 354, 163 

Pitkin M, Woan G., 2007, Phys. Rev. D, 76, 042006 

Prix R., Krishnan, B., 2009, Classical Quantum Gravity, 
26, 204013 

Punturo M. et al, 2010, Classical Quantum Gravity, 27, 
084007 

Santostasi, G., 2006, APS Meeting Abstracts, D1095 

Seto, N., 2005, Phys. Rev. D, 71, 123002 

Smits R., Tingay S. J., Wex N., Kramer M., Stappers B., 

2011, A&A, 528, A108 
Ushomirsky G., Cutler, C, Bildsten, L., 2000, MNRAS, 

319, 902 

The Virgo Collaboration, 2009, Adv Virgo baseline design 
Virgo Internal Note VIR-0027A-09 

Watts A. L., Krishnan B., Bildsten L., Schutz, B. P., 2008, 
MNRAS, 389, 839 

Will, C. M., 2006, Living Rev. Relativ- 
ity, 9, 3, URL (cited on 12 May 2011) 
|http : / /www ■ livingreviews . org/lrr- 2006- 3 | 

Willke B. et al, 2006, Classical Quantum Gravity, 23, S207- 
S214 

Xu R. X., 2003, ApJ, 596, L59 



